Transmission properties of a single metallic slit: Prom the subwavelength regime to 

the geometrical-optics limit 
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In this work we explore the transmission properties of a single slit in a metallic screen. We analyze 
the dependence of these properties on both slit width and angle of incident radiation. We study 
in detail the crossover between the subwavelength regime and the geometrical-optics limit. In the 
subwavelength regime, resonant transmission linked to the excitation of waveguide resonances is 
analyzed. Linewidth of these resonances and their associated electric field intensities are controlled 
by just the width of the slit. More complex transmission spectra appear when the wavelength of 
light is comparable to the slit width. Rapid oscillations associated to the emergence of different 
propagating modes inside the slit are the main features appearing in this regime. 

PACS numbers: 78.66.Bz, 42.25.Bs, 42. 25. Fx 



I. INTRODUCTION 



limits. Finally, in Sec. V we summarize our results. 
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After the appearance of the experiment of Ebbesen 
et al. [l| showing extraordinary optical transmission 
through a 2D array of subwavelength holes perforated 
in a thick metallic film, there has been a renewed inter- 
est in analyzing the optical properties of its ID an alog , 
an array of subwavelength slits 0, H S H E 0, 13, II ■ 
More recently, resonant transmission properties of a sin- 
gle subwavelength slit in a thick metallic screen have 
been analyzed both theoretically 'll] and experimen- 
tally Ha. Il4| . However, all these works consider the 
subwavelength regime. To the best of our knowledge, 
the study of the crossover between this regime and the 
geometrical-optics limit is still lacking. 

In this work, we explore how the contribution to the 
transmission current of the different modes inside a sin- 
gle slit evolves as a function of the incident angle and the 
slit width. To get physical insight, we apply a formalism 
based on a modal expansion of the electromagnetic fields 
already proposed in an equivalent problem 15, As 
we are interested in overall transmission behavior rather 
than in precise numerical comparison with experiments, 
we assume infinite conductivity for the metallic regions. 
Note, however, this approximation provides results of at 
least semi-quantitative value for good metals, such as 
silver or gold, if the dimensions of the structure are sev- 
eral times larger than the skin depthTo^. A bonus of 
using perfect metal boundary conditions is that all theo- 
retical results are directly exportable to other frequency 
regimes, simply by scaling appropriately the geometrical 
parameters defining the structure. 

This paper is organized as follows. In Sec. II we de- 
scribe the theoretical formalism. This approach is ap- 
plied to the subwavelength regime in Sec. Ill, where 
a comparison with available experimental data is also 
given. In addition. Sec. Ill presents two applications 
of a single subwavelength slit based on its characteristic 
resonant behavior. In Sec. IV we analyze the crossover 
between the subwavelength and the geometrical-optics 



II. THEORETICAL FRAMEWORK 

In this work we are interested in the scattering process 
between a plane wave and the structure sketched in Fig. 
1, i.e., a single metallic slit in a thick metal screen. Due 
to the translational invariance of the system along the y 
axis, it is possible to write the vectorial Maxwell equa- 
tions as two decoupled differential equations, one of them 
corresponding to s-polarization (electric field, E, parallel 
to y axis) and another corresponding to p-polarization 
(magnetic field, H, parallel to y axis). The case of s- 
polarization has been throughly analyzed in Here, 
we concentrate on p-polarized light in which resonant 
transmission properties have been reported. For this po- 
larization, the response of the structure can be solved by 
finding the solutions of the following equation 
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+ ^Hy{v) = 



(1) 



where we have assumed a harmonic time dependence of 
angular frequency uj. Hy{r) is the y-component of the 
magnetic field and r = {x,z). e(r) defines the spatial 
distribution of the dielectric constant. 

Once we know the field Hy (r) , the spatial components 
of the electric field can be obtained from 
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i dHy 
koe{r) dz 
i dHy 
fcoe(r) dx 



and. 



(2) 
(3) 



To solve Eq. (1) it is convenient to divide the space 
into three different regions (labelled as I, II and III in Fig. 
1). The magnetic field Hy can be expressed in regions I 
and III as a superposition of propagating and evanescent 
plane waves as 
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FIG. 1: Schematic plot of the structure under study. Dark 
grey corresponds to metalhc regions while light gray corre- 
sponds to a region defined by a homogeneous dielectric con- 
stant e. Both the reference system and the geometrical pa- 
rameters that define this system are plotted. 



H^ir) = Hi'"'{r)+ / dk p{k) exp[i{kx - k^z)lA) 
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where kz = \Jk'^ — k^, being fco the modulus of the 
wavevector of the incident plane wave. The coefficients 
p{k) and r(fc) are the transmission and reflection am- 
phtudes, respectively. The incident field is given by 
iJ*"c(r) = exp[i{koxX + ko^z)], being kox and koz the 
projections of the incident wavevector on the x and z 
axis, respectively (see Fig. 1). 

Inside the slit, the solution can be expanded as a sum 
of the eigenmodes of a waveguide 



^y^W = Yl eMiPmZ) + Bm{-il3mZ)](t>m{x) (6) 

m— 

where /?,„ = [fegC — (m7r/a)^]3 and the function 
4>,n{x) = {2/\^)coa[rmr/a{x + a/2)] corresponds to 
the solution of Eq. (1) within the slit, assuming perfect 
metal boundary conditions at a; = ±a/2. A homoge- 
neous material without absorption and characterized by 
a dielectric constant e, is considered inside the slit (light 
gray area in Fig. 1). By allowing e to vary we plan to 
check the reliability of a single slit as a fluidics detector. 

By applying the corresponding matching conditions to 
both Hy and at z = and z = W and proyecting 
the resulting equations on waveguide modes and plane 
waves, we obtain a system of linear equations with the 
set of coefficients {Am,Bm} as unknowns. We find that 
this problem is equivalent to that reported very recently 
in Refs. 15 and 16. In other words, the scattering prob- 
lem studied in this paper is equivalent to consider a fi- 
nite set of indentations with only one relevant mode in- 
side each indentation. Therefore, to obtain more physical 



understanding of the mechanisms that take part in this 
scattering problem, it is convenient to change the set of 
variables {Am, -Bm} by the new set {E^, E^} defined by 
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E° = ^ [Am exp {il3mW) - Bm exp {-il3mW)] (8) 
koe 

The set of equations to solve for the amplitudes 

{Em, Em} can be expressed as 

{fm + 9mm)Em+ ^ Qmrn' Em' + Qm^Z = -^m(9) 

{fm+9mm)Em+ 9mm" Em" + 9mEm = 0(10) 

m"^m 

We can interpret the equations (9) and (10) as follows. 
The rn-th mode has two associated field amplitudes, one 
of them (-E^) corresponds to the scattering events that 
take place at the input surface of the slit, while the an- 
other amplitude (E^) is related to the output surface. 
Each amplitude is determined as the result of different 
scattering processes between all modes. Now, we pass to 
analyze in detail the different terms appearing in the set 
of equations (9) and (10). 

Im is the overlap integral between the incident plane 
wave and m-th mode 



-^m = 2 / d:x.(j)m{x) exp{ikosin9 x) (11) 
Modes inside the slit are coupled through Qmm' 

dxdx'rm{x)cj>m'{x')G{±r,±'r') (12) 

The function G(r,r') = {TT/X)H^o\ko\r - r'|) is the 
two-dimensional Green's function corresponding to the 
vacuum, being Hq^\x) the Hankel function of the first 
kind. The unitary vectors x and x' correspond to x and 
x' axis, respectively. Notice that the propagator gmm' 
does not depend on the dielectric constant inside the slit 
but on the dielectric constant of the material outside, 
i.e., the coupling between modes takes places through 
the external medium. 

Amplitude of the m-th mode is inversely proportional 
to fm + 9mm, where fm is the associated admittance of 
that mode and Qmm is the self-interaction term. 
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expj-iPmW) + expjiPmW) 
(3m[exp{-i0mW) - exp{ipmW)] 



(13) 



where jm comes from the normalization factors of the 
wavefield inside the slit and it is equal to 1 if m = and 
1/2 if m ^ 0. 
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Amplitude of the m-th mode corresponding to the in- 
put surface is coupled to its counterpart in the output 
side through 



/3m(exp[i/?„Vr) - exp(-?/3„T4^)] 



(14) 



It can be proved that in the case of a non uniform slit, 
modes with different rn values are also coupled. 

It is possible to write the Fourier coefficients p{k) and 
T{k) as a function of {E^, E^} 



p{k) = 5{k - fcoc.) -Y^^'ra Sm{k) and (15) 

m=0 

oo 

r(fc) = Y.^^^rn{k) (16) 
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where 



Sm{k) — 



(17) 

Therefore the field Hy in air regions can be obtained 
from 



Hlir) = ijO(r) + £(-i?^) G„(r) (18) 

m=0 

oo 

H^'^ir) = J2EZGm{r) (19) 

m=0 

where we define Hy{r) = 2cos{kozz) exp{ikoxx) and 
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dx' Ct>m{x') G{t,t') 



(20) 



with r=(x,z) and being r' = {x',0) for the reflection re- 
gion {z < 0) and r' = {x' , W) for the transmission region 
{z > W). 

Equations (18) and (19) have a clear physical interpre- 
tation: scattered fields outside the sample are the result 
of the sum of those corresponding to a flat dielectric- 
metal interface [ff°(r)] and the contribution from each 
mode within the slit considered as an individual scat- 
terer. 

We can define the transmission current in the propa- 
gation direction as 



/oc 
dxJR{[E*(r) xH(r)]J (21) 
-oo 

If we normalize the transmission of the system (T) to 
the flux of energy that impinges on the slit we can write 
for region III 
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FIG. 2: Transmittance of a metallic single slit as a function of 
the wavelength of the incident radiation, computed with the 
same W=28.2 mm and two different values of a, 25 /im (a) 
and 75 fj,m (b). Notice that the wavelength range considered 
is the same in the two panels. Normal incidence and e=l 
inside the slit is assumed. 
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(22) 



As no absorption is present in the structure, current 
is conserved and the transmittance can be then alterna- 
tively and more easily evaluated inside the slit (region 
11) 



^=^^(E<*^™5-) (23) 



III. TRANSMISSION IN THE 
SUBWAVELENGTH REGIME 

First of all, we consider the transmission properties in 
the subwavelength limit, i.e., when X/a ^ 1. In order to 
compare our numerical results with available experimen- 
tal data, we begin by asumming W = 28.2 mm and the 
same range of wavelengths considered in Ref. 12. 

Figure 2 shows the transmittance in two diff'erent 
ranges of A/a within the subwavelength regime [in Fig. 
2(a) we take a=25 /xm while in Fig. 2(b) we consider 
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a=75 fim\. In this figure we assumed normal incidence 
and e—1 within the slit. In the subwavelength regime, 
only the propagating m — mode has a non-negligible 
contribution to the transmittance. Then, for this par- 
ticular case, the set of equations (9)-(10) converts into a 
system of two linear equations with a solution for Eq 
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{fo + gooy 



9f 



(24) 



It is then clear that the resonant transmission peaks 
observed in Fig. 2 are associated with the zeroes of the 
denominator of Eq. (24). In the strict limit a/X 0, goo 
goes to zero [as a/XLn{a/X)] and resonant transmission 
appears at the condition fo = ±3o, which straightfor- 
ward algebra shows equivalent to the usual Fabry-Perot 
resonant condition sin{koy/£W) = 0. When a/A is small 
but not zero, the transmission peaks are shifted and 
broadened due to the non-negligible contribution of the 
real and imaginary part of goo, respectively. Therefore, 
we can say that a subwavelength metallic slit behaves 
basically as a Fabry-Perot interferometer made of high 
reflective plates. Following with this analogy, it can be 
shown that when the slit width increases, the modulus 
of the reflection coefficients at the edges of the slit de- 
creases, and therefore the efficiency of the system as a 
interferometer becomes poorer. 

The numerical results shown in Fig. 2(b) are in agree- 
ment with experimental data published by Yang and 
S ambles [T^. However, it is worth noticing that for this 
case the expression obtained in Ref. 11, where only the 
leading terms in a/A are considered, overestimates (about 
10%) the shift of resonant peaks from the usual Fabry- 
Perot resonant condition, sm{kQy/eW) = 0, showing that 
the asymptotic a/A — *■ limit is reached very slowly. 

Another interesting feature observed in the subwave- 
length regime is that transmittance does not depend on 
the angle of incidence. This can be easily understood by 
looking at Eq.(ll). In this expression, the only quantity 
that depends on 9 is /q 



Io{0) = 2 sine 



(25) 



where sinc{x) = sin(x)/x. If a/A is small enough the 
dependence of this coupling on the incident angle is neg- 
ligible as sine « 1, independent on 0. 

Notice that, as it was mentioned before, we are nor- 
malizing the transmittance to the energy that reaches 
the slit. Therefore, since the total transmission in the 
resonant peaks for a=25/im is about three times that of 
a=75 /im, we can deduce that the energy flux collected 
by the slit at resonances is approximately the same for 
any value of a. This extraordinary property also implies 
that the electric field intensity inside the slit satisfies a 
1 /a dependence, and therefore, huge enhancements of the 
electric field intensity can be achieved just by reducing 
the slit width. However, as a leaky standing wave is ex- 
cited within the slit, only a periodic set of regions will 
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FIG. 3: Contour plot of the electric field intensity of a single 
metallic slit defined by a=40 fim, W=l mm. Incident radia- 
tion with A=409 fj,m and 9—0 is assumed. Red (blue) means 
maximum (minimum) field intensity. 



undergo the enhancement of the electric field intensity, 
as is shown in the contour plot of Fig. 3, where we have 
considered a representative example with a=40 ^m and 
W~l mm. Therefore, introducing a suitable non-linear 
medium inside the slit would produce a periodic mod- 
ulation of the refractive index in the direction of prop- 
agation. In other words, an one-dimensional photonic 
crystal is formed by the electromagnetic radiation at the 
resonant frequencies. 

It is also worth commenting that, as illustrated in 
Fig. 2, the linewidth of the transmission resonances de- 
creases when the slit width is reduced. This property 
allows us to measure the dielectric constant of small 
amounts of material filling the slit through the varia- 
tion of the resonant wavelengths with e, as has been sug- 
gested by Yang and Samblesfisl H^. In this paper, we 
study the accuracy of this kind of structure to work as a 
high-sensitivity microfluidics detector. Within this con- 
text, very recently, changes in transmission properties of 
photonic crystals have been proposed as a way to iden- 
tify certain very small quantities of fluids 1^. Here, we 
analyze a simple but functional alternative to this ap- 
proach. A small amount of a fluid can be guided to flow 
within the slit and just by recording the location of the 
transmission resonances, an estimation of the dielectric 
constant of the material can be obtained. To prove this 
proposal. Fig. 4a shows a contour plot of the transmit- 
tance as a function of the wavelength and the dielectric 
constant inside the slit. The efficiency of this system is 
shown in of Fig. 4b where the same two fluids consid- 
ered in Ref. 18 are assumed to flU the slit, i.e., e=1.8978 
(isopropanol) and e=2.2506 (xylene). As can be seen 
from this figure, the two peaks are well separated and. 
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FIG. 4: (a) Contour plot of the transmittance as a function of 
the wavelength and the dielectric constant inside the slit. Ge- 
ometrical parameters as in Fig. 2(b) are assumed, (b) Com- 
parison of the transmittance obtained with e=1.8978 (black 
line) and e=2.2506 (red line). 



in principle, with a slit with the geometrical parameters 
we have used in this simulation, this device could distin- 
guish whether isopropanol or xylene are filling the slit. If 
we were interested in characterizing two fluids with more 
similar dielectric constants, we could obtain even higher 
resolution in the transmission resonances by just reduc- 
ing the slit width. As commented previously, this device 
also is expected to work in the optical regime allowing 
then to characterize nanometric amounts of fluids. 



IV. CROSSOVER BETWEEN THE 
SUBWAVELENGTH REGIME AND THE 
GEOMETRICAL-OPTICS LIMIT 

In the previous sections we have found no dependence 
of the normalized transmittance with the angle of inci- 
dence in the extremely subwavelength limit. When the 
slit width increases, the dependence of the transmittance 
on incident angle becomes important. If a/ A is still small 
enough that only the first mode (m—O) contributes ap- 
preciably to the transmission, the coupling with the ex- 
ternal plane wave is governed by /q (6') [see Eq.(25)]. This 
is the case of the transmittance spectrum plotted in Fig. 
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FIG. 5: Transmittance computed with W=28.2 mm and (a) 
a=1.25 mm and (b) a=2.75 mm. Black, red, green and blue 
correspond to 9 =0°, 15°, 30°, and 45°, respectively. 

5(a), where we assume a=1.25 mm. As shown in this fig- 
ure, transmittance is maximum for 9 — 0° and decreases 
very slowly when 9 is increased. 

By increasing a, several eigenmodes with m > 
(both evanescent and propagating) start having a non- 
negligible contribution to the transmittance. The cou- 
pling of these m > modes with the external radiation 
and between them is given by /,„ and gmm' , respectively. 
An example of this behavior is shown in Fig. 5(b), where 
a— 2.75 mm. As can be seen in this figure, when a slit 
mode changes its character from evanescent to propagat- 
ing (dashed line in Fig. 5(b) shows the wavelength where 
TO = 1 becomes propagating, A = 2a), resonant transmis- 
sion properties change dramatically. Interestingly, this 
change is not observed for 9 — 0°. This is due to the 
fact that = 0°) = for TO > and gom — for 

TO odd. Then for this particular case of 9 = 0° and for 
this range of wavelengths, transmission is governed by 
mode TO = 0. Similar abrupt changes, at wavelengths 
where eigenmodes with higher values of to become prop- 
agating, have been also theoretically predicted for the 
transmission spectrum of light impinging on a circular 
aperture made in a perfect metal screen 19, 20] . 

Now, let us focus in the interval 0.1 < A/a < 0.275. 
Solid black line in Fig. 6 shows the transmittance in this 
interval of wavelengths for a=2.75 mm. For the sake of 
simplicity, in this section we analyze only normal inci- 
dence. If 9 7^0° similar physical mechanisms are present. 
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FIG. 6: Main: transmittance computed with W=28.2 mm 
and a=2.75 mm. Dashed lines show the wavelengths where 
the corresponding mode (see labels on top) becomes propagat- 
ing. Blue line corresponds to consider only the fundamental 
mode. Inbox: Comparison of the transmittance calculated 
taken into account all the modes (black line) and only m=0 
and m=8 (red line). Vertical green lines show the wavelengths 
where the Fabry-Perot condition is fulfilled for m=8 (see text). 

Dashed lines show the wavelengths where the correspond- 
ing mode changes it character from evanescent to prop- 
agating. First of all, notice that transmission peaks are 
associated to the emergence of propagating even modes 
and not to odd ones. This can be understood by real- 
izing that, as commented before, for normal incidence, 
odd modes arc uncoupled both to the external radiation 
and to the m=0 mode. Even modes, although uncou- 
pled to the normal incident radiation, have an indirect 
illumination through their coupling with m—0 mode. In 
addition, this makes that the fundamental mode approx- 
imately gives the overall transmission envelope (see blue 
line in Fig. 6). 

We have found that near the wavelength where the 
m-th mode (being m a even integer) changes its charac- 
ter from evanescent to propagating, i.e., at A/a = 2/m, 
transmission features are mainly governed by the con- 
tribution of only the mode m=0 and the m-th mode. 
In order to prove quantitatively this property, inbox of 
Fig. 6 shows the transmission in the neighborhood of 
the wavelength where m=8 becomes a propagating mode 



(A/a=0.25). Black line in this figure shows the result con- 
sidering all the modes within the slit, while red line shows 
the same calculation but now including only modes with 
771=0 and m=8 in Eqs. (9) and (10). We have found that 
771=0 gives the background of the transmittance while 
777=8 produces approximately the oscillatory behavior. 
These rapid oscillations of the spectrum can be explained 
by analyzing the z-component of the wavevector (fc^) as- 
sociated to the corresponding eigenmode when this mode 
changes from evanescent to propagating. Within the 
range of wavelengths in which this change occurs, /c^ — > 
and then Fabry-Perot condition sin{kz^/eW) = is ful- 
filled for many wavelengths inside this interval, resulting 
in transmission peaks with a very small wavelength sep- 
aration between them (see vertical green lines in inbox of 
Fig. 6). Finally, from Fig. 6 we can see how the ampli- 
tude of the transmittance oscillations decreases and tend 
to 1 when A/a 0, as it is expected from geometrical 
optics. 



V. SUMMARY 

In this work, a formalism for analyzing the transmis- 
sion properties of a single metallic slit under general con- 
ditions of angle of incidence and slit width has been pre- 
sented. 

First of all, we have carried out a comprehensive study 
of the subwavelength regime. We have discussed the ap- 
pearance of transmission resonances associated to the ex- 
citation of slit waveguide modes. We have also discussed 
two possible devices based on the huge enhancement of 
the EM-fields inside the slit that accompanies the ex- 
traordinary transmission phenomenon: a non-linear de- 
vice and a micro or nano-fluidics detector. 

In addition, the behavior of the transmittance when 
the subwavelength limit is abandoned has been also an- 
alyzed. We have studied in detail how the geometrical 
optics limit (T = 1) is recovered. In this regime, rapid 
oscillations in the transmission spectrum appear near the 
wavelength where a mode becomes propagating. 
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